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CHAPTER I 


INTRODUCTION 


In recent years communication system designers have begun 
to use the electromagnetic spectrum above 10 GHz in an effort to 
meet ever increasing bandwidth requirements. Overcrowding of the 
lower frequency bands have forced designers to make use of pre- 
viously unused portions of the spectrum. Higher operating fre- 
quencies are also necessary to support extremely fast data rates 
which accompany recent advances in solid state technology. 

At frequencies above 10 GHz, however, the effects of the earth's 
atmosphere upon communication signals dramatically increase. At 
these frequencies, a communications link designer must not only 
have a thorough knowledge of communication theory and systems de- 
sign but must also have a working knowledge of electromagnetic 
wave propagation. The necessity for this additional qualification 
is evidenced by the literature. Until just recently, the effects 
of the clear atmosphere fno precipitation) on a communications 
signal were assumed to be negligible. However, it has been demon- 
strated^ that clear air propagation effects can severely hamper 
a communications link under certain conditions. In fact, for sig- 
nal scintillations whose minimum levels approach the system margin, 
the probability of an error in detection increases significantly. 

Atmospheric effects upon a communications link can influence 
the amplitude, phase, angle-of-arrival , and polarization of the 
incident signal. These effects are manifest by amplitude and phase 
fluctuations, gain degradation, and depolarization. The causes 
of these effects include molecular absorption, scattering, refrac- 
tive index inhomogeneities and atmospheric noise. These causes 
are not, in general, in one-to-one correspondence with the effects. 

Molecular absorption at centimeter and millimeter wavelengths 
is primarily due to the presence of oxygen and water vapor. The 
former has a magnetic dipole moment, and the latter has an electric 
dipole moment. 

Scattering and absorption in rain at these frequencies is 
pronounced because the dimensions of the hydrometers are on the 
order of or larger than the operating wavelength. 

Refractive index irregularities arise from spatial and tem- 
poral variations in meteorological variables such as temperature, 
pressure, and water vapor content. These parameters are also in- 
fluenced by wind velocity and direction. 
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In clear sky situations, atmospheric gas emission is the pri- 
mary contributor to noise in this part of the spectrum. 

The major purpose of studying the effects of the earth's at- 
mosphere on electromagnetic waves at centimeter and millimeter 
wavelengths is to provide sufficient information for a communi- 
cations link designer to most effectively use this portion of the 
spectrum. Increasing research in this area has introduced new 
applications and techniques. 

In radio astronomy, for example, the effects of atmospheric 
propagation are important for several reasons. Earth-based radio 
astronomers must know what happens to a signal of extraterrestrial 
origin as it passes through the earth's atmosphere. Radio tele- 
scopes based outside the earth's atmosphere (such as those associated 
with Skylab) can use the knowledge obtained from the earth's at- 
mosphere to study foreign planetary atmospheres and vice versa. 

Recently, remote sensing has developed into a practical method, 
and in some situations the only method, of measuring meteorological 
variables at a distance. In particular, the structure constant 
of the refractive index and wind velocity perpendicular to the 
signal path can be calculated from statistics of a received sig- 
nal. 


One phase of experimental research on atmospheric propagation 
effects has been aimed at gathering long-term statistical data 
of a received signal. Statistics such as mean attenuation, fade 
distributions, and variance have been presented and analyzed in 
considerable detail. However, most of the analysis is performed 
in the time domain. The few papers that analyze the power density 
spectrum of the received signal estimate a slope and compare it 
with the theoretical value. 

In this thesis, spectral characteristics of 2 and 30 GHz sig- 
nals received from the Applications Technology Satellite-6 (ATS- 
6) are presented. These data were taken under clear air condi- 
tions. Second-order spectral characteristics such as slope and 
break frequencies are calculated at elevation angles ranging from 
0.38 to 43.89 degrees. The dependence of the spectral character- 
istics on elevation angle is investigated. 



CHAPTER II 


THEORETICAL BACKGROUND 


The earth's atmosphere is a random medium; therefore, wave 
propagation through the earth's atmosphere is a random process. 

As a result, literature dealing with this subject involves sta- 
tistics of atmospheric parameters. 

Centimeter and millimeter wave propagation is primarily af- 
fected by the troposphere. The ionosphere has little effect upon 
an electromagnetic wave at this frequency because the operating 
frequency is well above the plasma frequency. 

The first successful theory of microwave propagation through 
atmospheric turbulence was that of Booker and Gordon'. Difficul- 
ties in their theory gave rise to the^utilization of the statis- 
tical work of Kolmogoroff and Obukoff^*^*'^. 

Modern theories of atmospheric turbulence have physical foun- 
dations in the dynamics of fluids. The basic picture is that of 
a blob of fluid, a turb, which is put into rotating motion by tur- 
bulent eddies. The turb dissipates its energy either by perturb- 
ing surrounding turbs, which creates smaller eddies, or by gen- 
erating heat from frictional interaction with surrounding turbs. 

Large volumes of the fluid containing large eddies are as- 
sociated with significant kinetic energy whereas viscous frictional 
forces are thought to be smalr. Thus, the spectrum of atmospheric 
turbulence, $(u), as a function of spatial wavenumber, u, is di- 
vided into three sections: (1) the energy subrange, (2) the in- 
ertial subrange, and (3) the dissipation subrange' (see Figure 
1). The inner scale of turbulence, £j, is associated with the 
smallest eddy size, and the outer scale, £q, is associated with 
the largest eddy size. The inner and outer scales of turbulence 
are related to the spatial wave numbers, u. and u , as u. = 
and u^ 5 The inner scale of turbulence is on the order 

of 1 mm^. The outer scale of turbulence for an earth-space link 
can be taken as the effective height of the atmosphere since tur- 
bulent eddies larger than the effective height cannot exist. A 
value for the effective atmospheric height has been obtained by^ 
Devasirvatham and Hodge^ as fi km and will be used throughout this 
paper. 

The shape of the turbulence spectrum is different in each 
of the three subranges. Most papers dealing with this subject 
assume that form of the turbulencegSpectrum which is applicable 
in the intertial subrange given by^ 
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$(u) = 0.033 cl u'"' (1) 

2 

where C is the turbulence structure constant and m = -11/3 for 
the Kolmogoroff spectrum. 


$ (u) 



The validity of this assumption is based upon the size of 
the first Fresnel zone with respect to the inner and outer scales 
of turbulence, IaS., upper and lower bounds are found such that 
«n/~XL « where X is the operating wavelength and L is 
the effective signal path length through the atmosphere. 

Since propagation through the atmosphere is a random process, 
a signal received from a satellite is characterized in terms of 
covariance functions of the amplitude and phase scintillations. 
Expressions for the spectral density of the amplitude and phase 
scintillations are obtained by taking the Fourier transform of 
the covariance functions. 

Analytical theories start with Maxwell's equations or the 
wave equation and introduce scattering and absorption character- 
istics to obtain differential equations from which the covariance 
functions can be calculated. The expressions obtained are quite 
difficult to evaluate which leads to various approximate methods. 
On the other hand, transport or radiative transfer theories deal 
directly with the transport of energy through a medium containing 
particles. The basic integrodifferential equation is known as 
the equation of transfer. 


4 



Among the various approximate methods used, the most common 
is the Rytov method. The Rytov method assumes that the sum of 
the perturbations on the wave are much less than the magnitude 
of the incident wave and that the dimensions of all refractive 
index^perturbations are large compared to a wavelength. Lee and 
Harp‘d have shown that it is possible to relax these conditions 
by considering both the effects of these conditions on the mathe- 
matical expressions and the practical limitations of the measure- 
ment. Ishimaru^^ states that the Rytov solution is valid only 
when the variance of the log-amplitude fluctuation is less than 
0.2 vT 0.5. 

A commonly used approximation in the development of received 
signal covariance functions is Taylor's "frozen- in" hypothesis. 
Temporal variations in the received signal can be caused by re- 
fractive index changes along the signal path due to advection, 
convection, turbulent motions, etc. In addition, the motion of 
the atmosphere itself results in a transformation of spatial vari- 
ations of refractive index into temporal variations. Taylor's^ 
hypothesis assumes that the temporal variations of the refractive 
index are only due to the average wind velocity perpendicular to 
the signal path and, thus, can be removed by a simple transformation 
of the spatial variation, i.e.,8 

n(r,t) = n(7-vt) (2) 

where n is the refractive index, r is the position vector, V is 
the perpendicular wind velocity, and t is the temporal parameter. 
Thus, the effect of the frozen-in hypothesis on the time-lagged 
(temporal) amplitude covariance function, B (t), can be expressed 
mathematically as 

= <x(r,t) x(r,t+x)> 

= <x(r,t) x(r-"vT,t)> 

Bx(t) = C^(Vt) (3) 

where C is the spatial amplitude covariance function, t is the 
time la§, x is the log amplitude fluctuation, and < > denotes en- 
semble average. 

A quite significant assumption used in analyzing the ampli- 
tude fluctuations of the received signal is that of stationarity. 

The assumption of stationarity implies that the measurement is 
independent of the time it was taken. The validity of this assump- 
tion is tested in Chapter IV. 
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Mathematical Formulation of the Spectrum 
of the Amplitude Fluctuations 

1 q 

Using the Rytov method, it has been shown that the received 
field can be expressed as 

u(r,t) = u^(r)exp[t|;(r,t)] = u^(7)exp[x(r, t)+js(r,t)] (4) 

where u^(r) is the unperturbed incident field and the real part 
of \};, X, is the log-amplitude fluctuation and the imaginary part 
of s, is the phase fluctuation. 

The temporal amplitude covariance function is defined as 

Bx(t) = <x(r,t) xfr,t+x)> (5) 

where < > denotes ensemble average. Physically, B (t) is a measure 
of the similarity of two observations of the log-a(?lpl itude fluctu- 
ations separated in time by t seconds. If Taylor's hypothesis 
is assumed valid, the temporal amplitude covariance function is 
identical to the spatial covariance function (see Equation (3)). 

9 

Tatarskii has shown that the spatial amplitude covariance 
function for a plane wave propagating in a medium with a purely 
real refractive index and a Kolmogoroff turbulence spectrum given 
in Equation (1) takes the form 


Bx(t) = C^(vt) = 2tt / F^(u,L)J^(uvx)udu (6) 

where 




and where JQ(d) is the Bessel function of the first kind, v is 
the perpendicular wind velocity, L is the effective signal path 
length through the atmosphere, and k is the wave number (k=2TrA). 

The spectrum of the loq-ampl itude fluctuations is defined 
in terms of B (x) as a Fourier transform relationship, i.e., 

X 


F^(u.L) = irk^L 


1 - 


„2| 


SI 




^ •or* 

w^(f) = J dx 

— oo 

and since B (x) and W (f) are even 

X X 


( 7 ) 
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oo 


( 8 ) 


W (f) = 2 / B (t) cos{ 27 rfx)dT 

o’' 


Substituting (6) into (8) yields 


W^(f) = f 0.033 ^ 


11/6 


X J 1 - 

b ^ 


s i n ( ) 


( z+Q ^ ) 


exp 




j- 

z " dz 


(9) 


where 


f. = 


° 7 27tXL 


z = 


2, 

u L 


^ =7 27tXL ^ 


D = 


5.91 \ L 

^o J k 


Note that k and X are the operating wave number and wavelength, 
respecti vely, and are not related (directly) to the fluctuation fre- 
quency, f. 

Assuming that D»1 and changing the variable of integration 
to X = Equation (9) becomes 


2 ,7/6 , 11/6 Q 


W^(f) = j 0.033 k''^^ L 


-8/3 


X J 


1 - 


sin Q^{x+1) 


^ ( x+1 ) 


(x+l)”^^'^^ x"^^^ dx 


( 10 ) 


For small fluctuation frequencies Q « I and Equation (10) can 
be approximated as^ 
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, n «l 


( 11 ) 


W (f) ^ 0.075 <x^> ^ 

where <x^> = 0.077 ci For large fluctuation frequen- 

cies Q » 1 and Equation (10) is approximately 


W„(f) ^ 0.57 </>^ 


-8/3 


W^(f) = 0.57 <x^> Q. »1 


( 12 ) 


Equations (11) and (12) intersect at a break frequency, f|^ 2 ’ 
given by 



2.14fo 


2.14V 

v/2ttXL 


(13) 


Equation (13) provides a means for sensing the average wind veloc- 
ity transverse to the signal path. 

Figure 2 shows a Bode plot of W (f) for a mediom character- 
ized by a real refractive index. Thus, in a medium with a purely 
real refractive index, the spectrum of the log-amplitude fluctya^ 
tions is approximately flat up to rolls off as f 


W^(f) 



Figure 2. Theoretical spectrum of log-amplitude fluctuations 
in a medium characterized by a real refractive index. 
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Ott^^ has considered the spectrum of the log-amplitude fluctu- 
ations in a medium with a complex refractive index, n=a+jb. The 
spatial covariance function in this case is given by 

C^(vt1 = IIyP L / u$-j-(u)jQ(uvT)du H4) 

where $y(u) = 0.033 cl and cl is the structure constant 

of the temperature, t 1 The spectrurti of the log-amplitude fluctu- 
ations in this case is found from Equation (7) for low f as 


W^(f) = 0.002598 Cj (15) 

“ 8/3 

and exhibits an f** ' slope. 

Comparing (11), (12), and (15) an intersection frequency be- 
tween scattering and absorption regimes is found to be 


f 




(16) 


Thus, the spectrum of the log-amplitude fluctuations can be 
separated into regions dominated by a scattering mechanism and 
an absorption mechanism. This modified spectrum is shown in Fig- 
ure 3. 

The validity of this mathematical model of the spectrum of 
the log-amplitude fluctuations will be investigated in Chapter V. 

Wx(f ) 



Figure 3. Spectrum of log-ampl itude fluctuations in a medium 
characterized by a complex refractive index. 
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CHAPTER III 


EXPERIMENT DESCRIPTION 


The data presented in this thesis were obtained from the Appli- 
cations Technology Satellite number 6 (ATS-6) during the time period 
extending from P8 August to P5 October, 1976. The purpose of the 
experiment was to investigate the propagation of millimeter waves 
through the atmosphere, specifically to determine the effect of 
elevation angle on various communication system parameters. 

The receiving terminal was capable of measuring four frequen- 
cies: 30 GHz, 20 GHz, 2.075 GHz, and 360 MHz. However, the 20 
GHz transmitter failed just prior to the experiment, and the 360 
MHz data was unusable because of radio interference from local 
radio and television stations. Although only two frequencies were 
usable, their separation provided an excellent basis for data com- 
parisons. 

The ATS-6, a geosynchronous satellite, is used in some fifteen 
different experiments. Prior to the experiment reported herein, 

ATS-6 was located approximately 35^ East longitude over Lake Tan- 
ganika for the Satellite Instructional Television Experiment (SITE) 
based in India. In August, 1976, the satellite was moved to a 
position over the U. S. (^^94° West). This repositioning of the 
satellite provided an excellent opportunity for measuring the ef- 
fect of elevation angle on propagation characteristics. 

The satellite was moved an average of one degree per day and 
approximately one hour of data was taken at a number of different 
elevation angles except for the very low angles where almost con- 
tinuous data were taken. The observable elevation angles ranged 
from -0.48 to 44 degrees; however, angles below 0.38 degrees were 
considered unusable because of the possibility of ground reflec- 
tions. 

The satellite's antenna for 30 GHz transmission is a 0.46 
m parabolic reflector with linear polarization. The 2.075 GHz 
antenna is a 9.1 m parabolic reflector with right-hand circular 
polarization. 

The receiving terminal is located at the Satellite Communications 
Facility of the Ohio State University's ElectroScience Laboratory 
(ESL), 1320 Kinnear Road, Columbus, Ohio (latitude: 40°00'10"N, 
longitude: 83°02'30"W). The facility along with the antennas used 
in the experiment are shown in Figure 4. 
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Figure 4. Satellite communications facility. 



The 30 GHz receiver is shown in block diagram form in Figure 
5. The antenna is a horn Cassegrainean-fed 4.6 m parabolic re- 
flector with linear polarization. The beamwidth of this antenna 
at 30 GHz is 0.2 degrees. The front end (noise figure = 18 dB) 
consists of a solid state mixer and a stabilized local oscillator 
which produces a first intermediate frequency of 1.05 GHz. The 
front end is followed by a tunnel diode amplifier, manually-con- 
trolled step attenuators, and a Martin-Marietta phase-locked loop 
(PLL) receiver (modified for increased dynamic range). The re- 
ceiver bandwidth is 55 Hz, and the signal-to-noise ratio averaged 
52 dB at the higher elevation angles. 

The 2.075 GHz receiver was built at the ElectroScience Labo- 
ratory, and the block diagram of it is shown in Figure 6. The 
antenna is a horn front-fed, 9.1 m parabolic reflector with linear 
polarization. The beamwidth at 2 GHz is 1.3 degrees. The front 
end (noise figure = 8.5 dB) consists of a transistor mixer and 
a stabilized klystron oscillator which produces a first IF of 30 
MHz. The level of the 30 MHz IF was manually controlled by a step 
attenuator. The signal was then fed into the IF portion of a Collins 
75S-3 receiver to produce a second IF of 455 KHz, The receiver 
bandwidth is 4.5 KHz, which made phase locking unnecessary. The 
455 KHz signal was fed into a square-law detector. The average 
signal-to-noise ratio was approximately 52 dB. 

The outputs of the receivers were recorded on strip charts 
and input to a digital data acquisition system. 

Table I gives the link calculations for an elevation of 44 
degrees. 


Table I 

Link calculations (elevation angle = 44°) 



2. 075 GHz 

WWz 

Transmitter Power (dBm) 


33.0 

Spacecraft antenna gain (dB) 

39.5 

39.0 

Spacecraft system loss (dB) 

- 0.4 

- 1.0 

Free space path loss (dB) 

-189.9 

-213.4 

Gas Losses (H^O and 02)(dB) 

- 0.02 

- 1.4 

Ground antenna gain (dB) 

40.0 

57.8 

Polarization loss (dB) 

- 3.0 

- 

Ground waveguide loss (dB) 

- 

- 1.5 

Signal level at front end (dBm) 

- 70.8 

- 89.0 

Receiver noise temperature (°K) 

3800 

18000 

Receiver bandwidth (Hz) 

4500 

55 

Receiver noise level (dBm) 

-126.3 

-138.7 

Signal-to-noise ratio (dB) 

55.5 

51.2 
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30 GHz RECEIVER BLOCK DIAGRAM 



MARTIN MARIETTA 
PLL RECEIVER 


Figure 5. 30 GHz receiver block diagram 




2 GHz RECEIVER BLOCK DIAGRAM 
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2045 MHz 

COUPLER 


ATTENUATOR 


TO DIGITAL 
DATA 

RECORDING 

SYSTEM 


20 dB 


t 


KLYSTRON 
OSCILLATOR 
SRL 7G 
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RECORDER 


2045 MHz 



f© - 12 N f pgf ^ if ref 
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Figure 6, ?. GHz receiver block diagram. 




The data acquisition system is shown in block diagram form 
in Figure 7. The output of the receivers ranged from 0 to 5 volts. 
These analog signal levels were multiplexed and fed into an 8 bit 
analog-to-digital converter (A/D). A controller for the A/D sampled 
the data at 10 samples per second and 200 samples per second. 

The 10 Hz rate was always used, while the 200 Hz rate was recorded 
for approximately five minutes at each elevation angle. 

The main controller for the data acquisition system was an 
HP 2115A computer. It took care of sampling, formating, recording 
on 7 track magnetic tape, and keeping a log of operations such 
as attenuator settings, recording times, etc. The 7 track tapes 
were later converted to 9 tracks so that the data could be analyzed 
using the ESL Datacraft 6024 Computer. 

The analog signal levels from each receiver and wind speeds 
were also recorded on strip charts. 
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Figure 7. Digital data acquisition system 



CHAPTER IV 


THE DATA 


As mentioned previously, two sampling rates were used: 10 
Hz and ;?00 Hz. However, the data taken at the 200 Hz rate were 
not used in the work that follows because spectral frequencies above 
5 Hz were not significant and, thus, were for the most part con- 
taminated by system noise {noise in the spectra is discussed in 
the last section of this chapter). 

Figures 8 through 12 show some typical data samples at ele- 
vation angles of 0.38°, 4.95°, 18.11°, 33.65 , and 43.44 . Each 
waveform shown consists of 2048 data points taken at the 10 Hz 
sampling rate (1 data record = 2048 data points). There is ob- 
viously a significant effect of elevation angle on the received 
signals. This effect has been studied in the time domain by 
Devasirvatham and Hodge‘S. In this study, the effect of elevation 
angle on the signal spectra is considered. 

A. Stationarity 

15 

If the signal is "wide sense stationary" its mean, m, is 
independent of time and its autocorrelation function, R(x), is 
dependent only on lag time, t, i.e., it is not dependent on abso- 
lute time. In analyzing a data record, the mean of all 2048 points 
was calculated; thus, the mean is forced to be independent of time 
for any given data record (compensation processing for a nonsta- 
tionary mean is discussed in the next section). 

21 

The autocorrelation function is defined in general as 


R(t) = lim 

J->-co 


1 

2T 


T 


r 

-T 


f (t)f (t+T)dt 


n.7) 


It is more convenient, however, to calculate an autocovariance coef- 
ficient which is related to the signal waveform, f(t), by 


p ( T ) - lim 

J->oo 


T 


; 

-T 


[f ( t ) -m] [f ( t+T ) -m] d t 


(18) 


where m is the mean and a the standard deviation. For a periodic, 
sampled waveform. Equation (18) can be approximated as 
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TIME (SEC) 



7H I I 1 ! I 1 
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Figure 9, Received signals at el. = 4.95^. 
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Figure 10. Received signals at el. = 18.11°. 
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Figure 11. Received signals at el. = 33.65°. 
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o 



o 



Figure 12. Received signals at el. = 43.44^. 
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i=I 


a 


H9) 


with increasing accuracy as N-^co. 

If the autocorrelation function is dependent only on lag time 
then the time interval that is to be examined for stationarity 
can be divided into sections and the similarity of the autocorrela- 
tion function, and similarly the autocovariance coefficient, will 
determine the relative stationarity of the waveform. 

The data samples shown in Figures 8-12 were examined for sta- 
tionarity using the method described above. An additional 2048 
samples immediately following the samples shown in Figures 8-12 
were used so that the signal stationarity could be examined over 
the length of one data record. The first 2048 samples will be 
referred to as the left half and the second 2048 samples will be 
referred to as the right half. The autocovariance coefficients 
were calculated for both halves at each of the five sample ele- 
vation angles and are shown in Figures 13-17. For a numerical 
comparison, the squared differences between the left and right 
halves were defined as 


Squared 

Difference 


y n 


.yPL(ti) - PR(t^)]' 


rdBl 


7 .^P^t.^ . P^ft.)] 


^ 20 ) 


where p. is the autocovariance coefficient function for the first 
^048 points and pp is the autocovariance coefficient for the second 
2048 points. 

The squared differences for each of the five samples are 
shown in Table II along with the signal variances. Note that 
the stationarity of the sampled waveform improves steadily with 
increasing elevation angle. The nonstationarity at low elevation 
angles is probably due to large signal variances caused by strong 
fluctuations; as mentioned in Chapter II, signal (log-amplitude) 
variances greater than 0.2 to 0,5 are considered to be caused 
by strong fluctuations. The signal variances are shown in Figure 
16 as a function of elevation angle^. The curves represent re- 
gressive fits of the data to a spherical-earth model for the sig- 
nal path length. Ishimaru's strong fluctuation criterion cor- 
rgsponds to -7 dB in Figure 16. Thus, only the 30 GHz data under 
4° elevation angle is in the strong fluctuation regime according 
to this criterion. 
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Figure 14. Correlation between two consecutive 
data records at el. = 4.95°. 
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Figure 15. Correlation between two gonsecuti 
data records at el. = 18.11 . 
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Figure 16. Correlation between two consecutive data 
records at el. = 33.65°. 
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Table II 

Autocovariance coefficient comparison 


2 GHz Receiver 

El. (deg) 

Mean signal 

amplitude variance (dB) 

Squared 

difference (dB) 

0.38 

4.95 

18.11 

33.65 

43.44 

-17.5 

-25.5 

-37.0 

-42.0 

-43.0 

-14.3 

-18.6 

-32.6 

-38.7 

-42.6 


30 GHz Receiver 


0.38 

- 5.5 

- 6.1 

4.95 

-15.0 

-10.7 

18.11 

-28,0 

-22.5 

33.65 

-33.0 

-30.3 

43.44 

-40.0 

-35.2 


All of the correlation coefficients tend toward zero with^ 
increasing lag time except for the right halves at el, = 33.65 
and 43.44®. Their behavior is probably caused by short periods 
of no bit changes in the sampled waveform. 

B. Processing Techniques 

A Fast Fourier Transform (FFT) technique was used to approxi- 
mate the power density spectrum of the sampled signals fthe FFT 
actually produces a periodigram) . The periodigram was produced 
in two steps: (1) the FFT was used to obtain the Fourier Trans- 
form and its conjugate, F(w) and F*(o)), from the time-domain wave- 
form, f(t)'^; (2) then the conjugate pairs were multiplied and 
normalized by the d.c. component of the power. (The spectral 
components were normalized to the power in the d.c. component 
rather than the mean power so that the dependence of the spectral 
slope on variance was not removed.) The result, S(o)), was a per- 
iodigram which approximates the power density spectrum with in- 
creasing accuracy as the sampling rate increases. This two-step 
process is shown schematically in Figure 19. 


tNote that f(t) is in units of volts; thus, f(t) represents the 
amplitude fluctuation and not the log-amplitude fluctuations. 
However, the amplitude fluctuations and the log-amplitude fluctua- 
tions are similar if the fluctuations are small. 
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The highest spectral component that can be produced is the 
folding frequency, f^, which is half the sampling frequency, f^, 
i.e., f^ = f^/2. 

The lowest component is determined by the number of data points, 
N, and is found from the relation f = (f /2)/(N/?) = f^/N. For 
example, a period! gram produced from data sampled at a rate of 
10 Hz which uses 2048 points will contain spectral components at 
f = nf /N = 0.00488n Hz, where n=l ,2,3, ... ,1024. Periodigrams 
produced from the typical samples are shown in Figures 20-24. 

Two time-domain processing techniques were used to make the 
period! gram produced by the FFT more accurately approximate the 
true power density spectrum of the received signal: (1) slope com- 
pensation and (2) cosine tapering. The purpose of slope compen- 
sation is to remove sloping trends of the mean signal level. The 
effect of slope compensation on the periodigrams is that of re- 
ducing aliasing and Gibb's phenomena. Cosine tapering is a win- 
dowing technique the purpose of which is to reduce sidelobe levels 
in the spectrum. 

Slope Compensation Algorithm 

Sloping trends in the mean signal level are compensated by 
first considering M points fM^lOO for results that follow) on each 
end of the sampled waveform, see Figure 2S. Averages of these 
end samples are calculated, m^ and m^, and a straight line is drawn 
through these averages from P_. to P «. Each sample point is 
then level adjusted accord! ng^to the^offset level of the line through 
Pj^l and P^^. The slope compensation subroutine is given in Ap- 
pendix A. 
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Figure ?0. Per iodi gram of sampled waveform at e1 
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Figure 22. Peri odi gram of sampled waveform at el 
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Figure 23. Peri odi gram of sampled waveform at el 



2 GHz RECEIVER 



Figure 24. Periodigram of sampled waveform at el 






Figure 25. Slope compensation configuration. 


Cosine Taper Algorithm 

The window function, for the cosine taper is given 

+ cos T < t < T (21) 


by 


or in discrete form as 


1 


2 


1 + cos IT ^ 


f 5/\(vjp ( ^ -j ) » i “0 , 1 , 2 , . . . , M 


( 22 ) 


where t. is the time corresponding to the ith sample point and 
M is thJ number of points to be tapered on each end of the wave- 
form (the value of M was chosen as 20). The cosine taper subrou- 
tine is given in Appendix A. 

If f.-MoCt) is a pulse of width 2T then the Fourier transform 

of fQj(t) IS 


r f \ - irsinTo) 

r 1 CO / o 


0)(/-tVi 


(221 


-3 


Note that F^-p(o)) is oroportional to w**' as oh-oq 
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Figure 26 shows peri odi grams produced from each of these proc- 
essing techniques. In analyzing the data, both processing tech- 
niques were used to produce periodigrams; however, only 40 points 
(20 points on each end of a 2048-point data sample) were tapered. 

A third processing technique, filter compensation, was per- 
formed in the frequency domain. Referring to Figure 7 in Chapter 
III, note that the signals are low-pass-filtered before they are 
sampled. To compensate for the roll-off of this filter, each spec- 
tral component in a peri odi gram is increased by a value determined 
by the characteri Stic curve of the filter. Figure 27 shows the 
80 GHz filter characteristic curve and the effect of filter com- 
pensation on the spectral components. The effect of filter com- 
pensation on the 2 GHz spectrum is quite similar to that shown 
in Figure 27. 

C. Noi se 

It was desirable to determine the noise level with respect 
to the mean signal level so that the spectral shaoe could be ac- 
curately analyzed and interpreted. The noise present in the power 
spectrum can be classified into two types: receiver noise and bit 
noise. ^ Even though these types of noise are caused by different 
mechanisms, their effect on the power density spectrum cannot be 
separated. In this section, two methods are used to determine 
the noise level in the periodigrams. 

Method I. Consider the sampled waveform on a day when the 
signal was quiet (el. = 43.44^) as shown in Figure 28. A one-bit 
change on the A/D converter corresponds to 0.1,3 dB and is discern- 
able in the figure. Note that the A/D converter was toggling one 
bit above and below the mean signal level (0 dB). This random 
toggling could have been caused by receiver noise, bit noise, 
or the atmosphere. However, since a change of only one bit is 
observed, the receiver noise must be equal to or less than the 
bit noise. 

From the above argument, a waveform which is constant except 
for a single one-bit toggle should produce a spectrum whose maxi- 
mum component is a close approximation to the noise level in any 
given periodigram due to bit toggling. This level was calculated 
as -81.5 dB below the dc power level. 

Method II. Assume that the receiver noise is flat uo to 5 
Hz and drops to zero above 5 Hz. Then with a system margin of 
^2 dB and spectral lines spaced 0.00488 Hz apart, the noise level 
is calculated as 
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Figure 27. Low-pass filter characteristics and appropriate 

spectral compensation. 
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Figure 28. Sampled waveform on a "quiet" day. 


— ^ - 52 = -82.0 dB (24) 

A peri odi gram of the sampled waveform shown in Figure 28 is 
shown in Figure 29. A linear regression line drawn through the 
spectral components intersects the ordinate at -76 dB. 

Considering the noise levels given above, a threshold level 
at which the noise "corrupts" the spectrum was first set at -80 
dB. However, after examining the results of slope and break fre- 
quency calculations (discussed in the next chapter), this "noise 
cutoff level" was raised to -70 dB below the dc power level, 
which corresponds to the peaks of the spectral components in Fig- 
ure 29. 


noise level = 10 
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Figure 29. Periodigram of “quiet" sampled waveform. 
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CHAPTER V 


RESULTS 


In this chapter the spectral characteristics of the experi- 
mental data are investigated. A comparison is made between the 
experimental results and the theoretical results discussed in Chap- 
ter II. The spectra are characterized by either one or two slopes 
which are obtained by a linear regressive (least squares) fit. 
Probability density functions are given for slopes and, in the 
case of two-slope fits, break frequencies. The dependence of the 
spectral shapes and break frequencies on elevation angle Ipath 
length) is also investigated. 

A. Approach 

The basic approach was to process the time-domain digital 
waveform using the techniques described in Chapter IV to obtain 
a reasonably accurate power density spectrum. 

The next step was to fit the spectra to a shape similar to 
that shown in Figure 3. Note that the lower break frequency in 
Figure 3 is dependent on the rate of change of the real and imag- 
inary parts of the refractive index with respect to temperature, 
9a/BT and 9b/9T respecti vely. The values for these quantities 
are calculated in Appendix B using oxygen and water vapor as the 
absorptive gases. These values were used in Equation (16) to com- 
pute the lower break frequency, f . , . The results of these calcu- 
lations for 50% humidity at an elevation angle of 20° are given 
in Table III. The values of f. , given in Table III are. several 
orders of magnitude lower than°those calculated by Ott"^^ using 
oxygen as the only absorptive gas; the reason for this discrepancy 
is unknown at this time. 


Table III ^ 

Lower Break Frequency (el. = 20^) 


Frequency fGHz) 

aa/9T/9b/9T 

fbi (Hz) 

2.075 

30.000 

R.n? X 10^ 
fi.8o X in' 

4.71 X 10'^^ 
2. on X in' 


The above calculation indicates that the lower break frequency 
should not be seen in the periodi grams. Therefore, a maximum of 
two slopes were fit to the spectra. 
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Two slopes were fit to the spectra by an iterative procedure 
based on the minimum of the root-mean-squared (r.m.s.) errors of 
each regression. The first fit, characterized by a slope, Al, 
and intercept, P., was obtained by setting the break frequency 
(the frequency at which the two slopes intersect), f. , equal to 
the eleventh harmonic*, performing a linear regression on the first 
eleven harmonics and then performing a second linear regression, 
characterized by a slope, A2, constrained about the point at which 
the first regression curve intersected a vertical line drawn through 
the break frequency, see Figure 30. 


S(f ) 



Figure 30. Geometry of two- si ope fit. 


Each successive fit was obtained by incrementing the harmonic 
number and fitting two more slopes. After obtaining the slopes 
for the maximum harmonic number, r.m.s. errors at each break fre- 
quency were compared. The slopes obtained with the minimum r.m.s. 
error were considered the best fit, and the harmonic frequency 
at which the best fit was found was taken as the break frequency. 


*The eleventh harmonic (O.OF37 Hz) was chosen as the first usable 
break frequency to ensure that at least ten points were used for 
the first regression. 
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By visually examining several spectra, the maximum harmonic 
number for the break frequency search was set at 250, which cor- 
responds to break frequency harmonic was not 

allowed to come^closer than ten harmonics within the harmonic cor- 
responding to the noise cutoff frequency. The break frequency 
search algorithm along with the noise cutoff algorithm is given 
in Appendix A. 

After obtaining two slopes and a break frequency, a noise 
cutoff level adjustment was made to prevent the slopes from being 
contaminated by the noise in the spectrum discussed in Chapter 
IV. From the results presented in the last section of Chapter 
IV, the level at which the noise corrupted the spectrum was -70 
dB. Therefore, the levels at each point on the regression lines 
for the spectral slopes were compared with the noise cutoff level. 

The frequency at which the calculated slope lines intersected the 
noise cutoff level line (point S in Figure 30) was termed the noise 
cutoff frequency, f . Once f was determined, the entire break 
frequency search procedure was repeated; however, only the harmonic 
numbers below that corresponding to f^ were used in the regressions. 

Figures 31 through 35 show two-slope fits as calculated from 
typical samples. In Figures 32a, 33a, and 34a it is obvious that 
the two slopes do not accurately represent the spectrum. The prob- 
lem is most probably in the error minimization criterion and in 
allowing the break frequency to come too close to the noise cutoff 
frequency. 

Because of the problems with the break frequency search routine, 
the results were visually classified into four types: (1) A1 > A2, 

(?.) A1 < A2, (3) A1 A2, and U) other unacceptable cases due 
to difficulties in the computer algorithm. The type (3) cases 
were recalculated using a single linear regression curve fit (the 
single slope regression algorithm is given in Appendix A). Only 
type Q) cases were used for break frequency statistics. Slope 
statistics were found from the second slope of type (1) cases, 
and the first slope of the type (?) and (3) cases. The second 
slopes of the type (2) cases were assumed to be caused by noise 
in the spectrum and, thus, were not used. All type (4) cases were 
considered unusable. The statistics of the data sets are given 
in Table IV. 
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(b) 

Figure 31. Two-slope fit at el. = 0.38° 
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Figure 32 . Two-slope fit at el. = 4.95^. 
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Table IV 

Results of Data Classification 


Type 

No. of Cases 
at 2 GHz 

No. of cases 
at 30 GHz 




1 

142 

238 

2 

IIP 

53 

3 

112 

117 

4 

98 

62 

Total 

471 

470 


B. Probability Density Functions 

In Chapter II it was, shown that the spectral slopes should 
roll off according to f”°' , which corresponds to -26.7 dB/decade 
on the peri odi grams. In this section, the distribution of the 
slopes and break frequencies are described by probability density 
functions. Figures 36 and 37 show the slope probability density 
functions. The mean of the 2 GHz data is -23.5 dB/decade and 
-21.8 dB/decade for the 30 GHz data. The distributions about these 
means are quite different. The 2 GHz data has a very high peak 
at -25 dB/decade and has a wide range. The 30 GHz data does not 
contain a single high peak as seen in the 2 GHz data; in fact, 
the 30 GHz data exhibits a double peak behavior. This type of 
behavior might be the result of bit noise contamination, i.e., 
the density to the left probably corresponds to turbulence- 
induced spectra while the density to the right might be the result 
of spectra contaminated by bit noise. 

There are several possible explanations for the average value 
of slope to be less steep than the theoretical value. Random bit 
toggling would cause the spectral slope to approach -20 dB/decade 
or flatter (see Appendix C). Aliasing would tend to enhance the 
higher frequency components which would also make the slopes less 
steep. Causes of the variations in slope might include spatial 
and temporal inhomogeneous turbulence, changes in average wind 
direction, and nonstationarity of the process. 

The break frequency probability density functions are shown 
in Figures 38 and 39. Note that the bulk of the frequencies is 
at the lowest value that can be obtained from the break frequency 
search algorithm. This high density toward the lower values sug- 
gests that there might be a significant number of break frequen- 
cies just below the minimum value. Also, since the wind speed 
data have a fairly symmetric distribution, the break frequencies 
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Figure 36. Probability density function of spectral slopes 

for 2 GHz receiver. 
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Figure 37. Probability density function of spectral slopes 

for 30 GHz receiver. 
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Figure 38. Probability density function of spectral 
break frequency for 2 GHz receiver. 
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Figure 39. Probability density function of spectral 
break frequency for 30 GHz receiver. 
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should also have a fairly symmetric distribution. Thus, some of 
the break frequencies might represent second order break frequencies. 
For example, if the actual break frequency was below the minimum 
(0.05 Hz), then the break frequency search would attempt to find 
a break frequency higher than 0.05 Hz. 

The range of values of break frequencies is quite reasonable 
when compared with break frequencies corresponding to appropriate 
wind speed data. For example, a break frequency of 0.1 Hz for 
the 30 GHz data corresponds to a perpendicular wind soeed of 1.54 
m/sec at an elevation angle of 20°. For the 2 GHz data, a break 
frequency of 0.1 Hz corresponds to a perpendicular wind speed of 
2.38 m/sec. 

Since wind speed was recorded on a strip chart and wind speed 
and direction was also available from the weather bureau (about 
4 miles east of the receiving terminal) a comparison between meas- 
ured wind speeds perpendicular to the signal path and wind speeds 
calculated from the break frequency data was made. However, the 
correlation was quite low, particularly at the low elevation angles; 
the wind speeds calculated from break frequencies were mostly greater 
than those measured directly. A poor correlation was not surprising 
since the measured wind speeds were point measurements on the ground 
and the calculated wind speeds represent an average value along 
the signal path. For this reason, inhomogeneities in the atmos- 
phere could cause poor correlation between measured and calculated 
wind speeds, especially at the low elevation angles. 

C. Dependence of Spectral Characteristics 

on Elevation Angle 

According to the discussion in Chapter II, the spectral slopes 
should obey the f” ^ (-26,7 dB/decade) law. No dependence on 

elevation angle is found in the weak fluctuation theory. However, 
as mentioned earlier, the 30 GHz data at low elevation angles (below 
4°) does not satisfy the weak fluctuation criterion according to 
Ishimaru. The strong fluctuation theory predicts a spectral slope 
proportional to f or 0 dB/decade (see Appendix D). It is not 
clear at this point how the strong fluctuation theory interfaces 
with the weak fluctuation theory. 

Figures 40 and 41 show the dependence of spectral slope on 
elevation angle. The crosses represent maximum and minimum values 
and the circles are the average slopes, in dB/decade, at each ele- 
vation angle. The dashed line indicates the theoretical -26.7 
dB/decade slope, and the solid line represents the average value 
of all measured slopes. The effective path length scale is based 
on a spherical earth model with a 4/3 radius and a homogeneous 
atmosphere of 6 Km. The effective path lengths ranged from 8.6 
Km to 318 Km. 
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Figure 40. Dependence of spectral slope on elevation 
angle for 2 GHz receiver. 
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Figure 4t, Dependence of spectral slope on elevation 
angle for 30 GHz receiver. 
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The 2 GHz data has a quite large maximum-to-minimum spread, 
although the average is quite close (3.2 dB/decade above) to the 
theoretical value. The average points at the low elevation angles 
and at the high elevation angles are in good agreement with theory. 
Most average points between the elevation angle extremes are sev- 
eral dB/decade higher than the theoretical value. 

The 30 GHz slopes have a much narrower maximug-to-minimum 
spread. Data taken by Vogel, Straiton and Fannin‘S” is also shown 
in Figure 41. Both data sets seem to be in fairly good agreement, 
with the data taken by Vogel, et al. tending to have less steep 
slopes at the lower elevation angles. It is unknown whether the 
fact that the data obtained by Vogel, et al. was measured over 
water might have caused less steep slopes at low elevation angles. 

Note that there seems to be a trend toward steeper slopes 
from el. = 0.38° to el. = 5°. This trend is also apparent in the 
data taken by Vogel, et al. This trend might have something to 
do with strong fluctuations; however, the slopes do not approach 
0 dB/decade as predicted by the strong fluctuation theory. Aper- 
ture averaging might also be an influence of this trend since aper- 
ture averaging would reduce the small-scale fluctuations which 
would in turn cause steeper slopes. 

Immediately following 5^, the slopes show a sudden change 
toward less roll off. The data taken by Vogel, et al. also shows 
this sudden change, but it is not as pronounced as the O.S.U. 
data. This sudden change in values seem to suggest a transition 
region between two mechanisms; for example, between strong and 
weak fluctuations or between angle-of-arrival fluctuations and 
amplitude fluctuations. 

The slopes above 5° are mostly flatter than the theoretical 
value. The reason for this is again probably due to bit noise. 

The 2 GHz break frequencies as a function of elevation angle 
are shown in Figure 42. The spread is seldom more than 0.05 Hz, 
and, as pointed out in the p.d.f. 's, most of the break frequencies 
are close to the lowest possible frequency which could be obtained 
from the break frequency search algorithm. 

The 30 GHz break frequencies as a function of elevation angle 
are shown in Figure 43. The spread is much more than the 2 GHz 
data, sometimes as much as a factor of 20 greater than the 2 GHz 
data. The 30 GHz break frequencies are not as heavily concentrated 
about the lowest value as would be expected if the break frequen- 
cies are indeed related to the perpendicular wind speeds as given 
in Equation (13) (Chapter II). 
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Figure 42. Break frequency as a function of elevation 
angle for 2 GHz receiver. 
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D. Comparison of Spectral Characteristics 
between 2 and 30 GHz Data 


To gain a better insight as to the nature of the data, the 
spectral characteristics at 2 GHz and those at 30 GHz were compared 
visually from scatter plots. Figure 44 shows a scatter plot of 
the 30 GHz slope versus the 2 GHz slope. Theoretical ly, the spec- 
tral slopes should be -26.7 dB/decade at both frequencies. This 
scatter plot shows that the slopes are indeed heavily concentrated 
around the theoretical value. 

Figure 46 shows a scatter plot of the 30 GHz break frequency 
versus the 2 GHz break frequency. Theoretically, the ratio of 
the 30 GHz break frequency to the 2 GHz break frequency is 3.87, 
as indicated by the dashed line. The solid line represents a linear 
regression through the scatter points. The agreement with theory 
is quite close even though the lowest attainable break frequency 
is most likely too high. 
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Figure 44. Scatter diagram of 30 GHz spectral slope 
versus 2 GHz spectral slope. 
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Figure 45. Scatter diagram of 30 GHz break frequency 
versus 2 GHz break frequency. 
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CHAPTER VI 


SUMMARY 


The results presented herein characterize spectra of ? and 
30 GHz sianals received from the ATS“6 at elevation anqles ranqinq 
from 0.38^ to 44^. 

The spectra of the received signals are characterized by a 
two-slope linear regressive fit and a break frequency corresponding 
to the intersection of the two slope lines. Cases in which the 
two slopes were approximately equal were recalculated for a single 
slope fit. The results are presented statistically in the form 
of probability density functions. The dependence of the spectral 
characteristics on elevation angle is given. Scatter diagrams 
are used to compare the 2 and 30 GHz results. 

Average spectral characteristics are in general agreement 
with predictions based on turbulence theory; however, variations 
of the slopes can be quite large. In fact, temporal variations 
of spectral slopes are vecy similar to variations of the refractive 
index structure constant"^'. Possible causes for variations in 
slope include spatial and temporal inhomogeneous turbulence, changes 
in average wind direction, and nonstationarity of the process. 

A trend toward steeper slopes appears to be evident at low 
elevation angles (below 5°) in the 30 GHz data. The cause and 
certainty of this behavior is unknown at this point and should 
be carefully examined in future experiments. 

Because the time-domain waveforms are mostly low-order bit 
toggles at elevation angles greater than 25°, the spectral slopes 
and break frequencies above 25° are probably biased representa- 
tions of the spectral characteristics. 
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APPENDIX A 


COMPUTER ALGORITHMS 


the 

The following routines 
techniques discussed in 

are computer algorithms which 
this thesis. 

implement 


Technique 

Subroutines used 

Lines 

1. 

Slope compensation 

MEAN 

1-40 

P. 

Cosine taper 

COSTAP 

44-60 

P. 

Power density spectrum 

FFT 

64-106 

4. 

Break frequency search 

CORRUPT, BREAK 

111-18? 

5. 

Single slope regression ONESLP 

186-213 

6. 

Regressions 

LINREG,CLR 

217,298 
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It. c CORRECT FOR shifts OF MF/\N LfVel BY SAMPLING NeND POINTS 

I7 r FROM TTMe-1!OMaTN WaVFFORm 

10 AVP=0, 

19 MP0INT=?n4A-NEM0 

P-8 DO 10 Isl.NEND 

21 AV1=AV1+Y(I) 

K=:MPOlNT4l 

23 in AV? = AV;> + Y (K 1 

24 AVt=AVl/NEMI; 

25 AVP=AV2/NEND 

26 AVniF= ( AVa-AVl ) / ( 2fi4fl«MEND ) 

27 AVCOr\}=AVl-AVniF'MNEMIV2. ) 

28 00 20 I=l,pn4A 

29 ?0 Y< I )=Y C I )-AVcOM- ( AvijIf^t ) 

■50 A\/1 = ( AV1+AV2 ) /?» 

31. PflTURM 

32 C CAI CULATE MEAN OF SaMPI EO OaVFEORM aNU SUOTRaCT 

33 r FROM Y 

34 30 no 40 1=1.2040 

35 40 AVlsAvl+Y(l) 

36 AVl = AVl/2t)40 

57 00 50 1=1,2040 

56 50 Y<U=Y(I)-AV1 

59 return 

4(1 FND 

41 c 

42 r. 

43 r 

44 SUPROUTINE rOSTAPlY.M.N) 

4b r 

46 C COSINE taper 

47 C TIME-DOMAIN SAMPLED WaVEfUF^M 

46 r MsMUMBFR (»F SAMPLES T() Rr TAnFKfi) ON EACH ENii OF RECORD 

49 r m=MUMRFr of SAMPLES 

b (1 C 

bl n I MENS I ON Y(l> 

b2 A=PH1.) 

b3 I'O 10 T=J,M 

54 vIsN + t-T 

bb )/M 
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56 r=.5*(l.+C0S(R) ) 

57 Y(T)=Y(I>*C 

5(\ 10 Y«.)) = Y(J)*r 

59 nt:TURfj 

60 Ff\jn 

fi-i r 

r, 

6^ r 

6m SURROUTIME FFT<X»Y,S) 

65 c 

66 r PHnOUCF PfRiOOlGRAM FROM SAMPLER WAVEFORM 

67 C nSiNG FFT (FORT) 

60 C 

69 r X -- ARRAY OF HAHmonIC FreQUFMCIFS 

70 C (MULTIPLES OF TLF) 

71 r Y -- TIMt-l.iOMAlM SAMPLED WAVEFURm 

7'^. r .9 -- array of spectral CoMPFMENTS 

75 C 

7M complex B 

75 COMMOiM/ARAYS3/P(iiO‘F0) «SS(512) 

76 OIMENSIOM X« J02‘») «Y(2OH0) *S(10;;J4) 

77 COMMOU/SLll /fJFR£« t MRATF ♦ TLF 

70 C0MM0i\i/SET5/UnTTP^ match »MtAP » KAv.CUtOF « AVI t FTP 

7P nATA rjPMTS/204A/»NHF/102M/ *M/H/ 

60 TLF=FLOAT(NL'AtE) /float (NpNTS) 

01 r '^ET UP SINE Taple_ 

62 CALL F0RT(R,M,SS,0,IFFRR) 

05 no 5 I=1.NPNTS 

8m p ( I )=CMPLK ( Y ( I ) *0 . ) 

65 6 COMTINUF 

06 r FAST fourifR transform SaMPlfo waveform 

07 CAl.L FnRT(BtM,SSi~2,IFFRR) 

80 TF ( IFEE^E(.NE»0 )MlHR=Y( P( 15HF RRDH If'J FORT, 15.0) 

89 r CAI CULATE OC COMPONENT 

9U RA = REAL (0(1) )**2 

91 nCP=AVl*AVl+RA 

92 C calculate POWf.R SPECTRaL COMPONENTS 

95 C IN DR BELOW M( AN LEVEL 

9M no 10 T= 2 .nHF 

95 jsNPMTS+^-I 

96 L=:I-1 

97 SW = 2.*REAL(EMI »*B( J) ) 

90 p(L)=10,*AI 0(iJO(SW/OCP) 

99 XC=TLF*L 

ion X(L)=ALOG10<XC) 

101 10 continue 

102 r CORRECT spectral COMPONENTS FOR 

105 C LOW-PASS filter KULL-OPF 

lOM call rwcoro.s) 

105 RETURN 

106 ENH 

107 C 
106 r 
109 r 
Ho r 
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Ill 
11? r 
11^^ c 

114 c 

iis r 
li<-- r 
117 c 

lift c 

119 
12(1 
121 
122 
12.^ 

1?4 

12 fi 10 
126 
127 

120 ?0 
129 

1 JO 
IJl 
13;> r 
1J3 r 
134 r 

1JT1 
U6 r 
137 c 
150 r 
159 r 
14U r 

141 r 

142 r 

14? c 

144 r 

145 c 

146 r 

147 r 
14 m 

149 

150 

151 

152 
15? 

154 

155 

156 10 

157 

158 

159 

160 r 

161 
162 f 
16? 

164 r 
16b 


<;UrROUTII\IE CORUPTCcLlFtCnTOD 

WOISE r.UTOFF AD^USTMFNTI 

FIMO NOISF CuT< 1FR FRfc.OuPMCY ANj SEARCH 

FOR RREaK FREOHENCY ONLY UP TO tmAT FREOUEMCY 

CLIF -- NOISE CUTOFF FRfQUfMcY (HZ) 

CUTOF -- NOISE CUTOFF LF\/EL ( OH ) 
(:0HM0N/SET4/FB»AI,P1, A?»rR 
rOI^MON/ARAYSl/X( 1024) , Y ( ?o40 ) » S () 024 ) 
COMnON/ARAYSp/XRl 1024 ) »Yp(2n«ifl) ,Y0UT(l0?4) 
nATfl NMFM1/1023/ 

TLF=:10,/2040, 

no lo i=i,NnFMi 

TF(YOUT(I),LT,rUTUF)GO TO 20 

CLTFsS. 

return 

CLIF=I*TLF 

call BREAK(y,S»FB,Al«Pl»A2»FPi I ) 

return 

FNO 


SuPROuTTNE break ( X , YtFP, A1 ,P1 , Ap , ) 

PRFAK FREOUFMcY search Afll) TWU-SL OPf HT 

Y -- array oF harmonic FRF(?UEf,iC It-S 

Y ARRaY of spectral COC'ipOiMENTS 
FR -- break FREOUENLY 

A1 -- FIRST SLOPE (DR/nrC) 

8t -- power level intercept <DB) 

A2 -- SECOND slope (DFI/PEC) 

ER -- total RMS ERROR 

mend-- ending harmonic For sfahc)< 

dimension xa)tY(i) 

COMMON/SET1/NFREO*NRATE»TLF 

COMMDN/ARAYS2/F(in?4) ,YP(204O) «Y0UT(lDp4) 
data NMIN/10/,MMAX/250/ 

TLF=10. /204ft. 

MFR=NMIM 
IMTN=1 
1 = 0 
T = I4-1 
NFB=NFB+1 

noif=nenu-nfb 

IF (NDIE ,LT,10)G0 TO 15 

calculate first slope AMp POWER LEVEL INTERCEPT 
CALL LTMREG(X.Y»1»MEH, Al ,E^1 »fRRi , YOUT ) 

CAICULATE SfCONC) SLOPE 

CALL CLr ( X t Y »NFR»NI ND , A2 , E RP 2 • YODT ) 

SLIM RMS ERRORS EHOb ROTH RE OR F SS I OIjS 

r( J)=ERHI+KKH2 
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167 

16f> 15 

169 r 

170 

171 

172 
1 73 

174 20 

175 r 

176 

177 
176 
179 
IHO 
IBl 
182 

183 r 

184 r 

185 C 

18b 

187 C 

188 r 

189 C 

190 r 

191 r 

192 C 

193 r 

194 C 

195 r 

196 r 

197 r 

198 

199 

200 
201 C 
2t»2 

203 r 

204 

20 5 10 
206 

207 

208 20 

209 C 

210 C 

211 30 

212 

213 

214 C 

215 r 

216 r 

217 

218 r 

219 r 
22ii r 


TF(NFB.RE.,IMMflX )60 TO 15 

60 TO lO 

FRsLd) 

find minimum tRROR 
no 20 .1 = 2,1 

IF<F( J) .GT,FR)GU to 20 
FH=E( J) 

IM.TN=J 

rOMTINUE 

9 ECAI-CULF 0 RF. SLOPES FOR MINIMUM f RKOR 
MFP=WMIM+IMIN 

tai l LIMREGf X«Yd •MFP«A 1 ,B 1 iFRHi ,Y 0 UT) 
call CLR(x, Y,NFB,NEIMn,A 2 ,ERP?« YOUT) 

FB = NFE)*TLF 
rR=ERKl+ERR 2 
Return 
run 


subroutine ONFSLPly ,S,Al,Bl,rR,Ct IF,ruTOF ) 
SINGLE SLOPF REGRESSION 


X 

S 

Al 

Bl 

E.F? 

CLIP 

rUTOF 


ARRAY OF harmonic FPEOliENCiES 
ARRAY OF spectral COMPONENTS 
SLOPF (UB/OLC) 

POWER LLVFL intercept (OR) 
TqTAl RmS ERROR 
NOTSF cutoff FRfOUEmCy (HZ) 
MoISf cutoff level (UB) 


n I mens I ON X ( 1 ) t S( 1 ) , Y0UT(102ti ) 

OATA NHFMI/1023/ 

TLFs:10./204n. 

CALCULATE SLOPE 

fall LlNREG(X,Stlt]023.Al,B3 ,ER,YOUT) 

FIND NOISE CUTOFF FREQUENCY 
no 10 I=1 ,NHFm 1 
TF(Y0UT(I),LT,CUT0F)G0 To 20 
CLIF=5. 

GO TO 30 
CLIF=I»TLF 

RECALCULATE SLOPE USING SPECTRAi COMPONENTS 

OELOW noise cutoff FREQUENCY 

CALL LlNREG(X,S»ltI,Al,Ol ,ER,yOuT) 

Return 

FNO 


subroutine LINPEG(X»TfNLO»NMT»SLOPE,YlNT,tR,YO|jI ) 
I.INLAR ULGRFSSTON 
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2P1 C 
22 ? r 

223 c 

224 C 

225 C 

226 r 

227 r 
220 C 
229 c 
250 C 
231 
252 

233 

234 

255 

256 
23/ c 
23fi 

239 

240 

241 
24? 3 

243 

244 

245 
?4(^ 
247 
24B 

249 C 

250 

251 

252 

253 2 

254 

255 

256 

257 
256 c 
259 C 
2fe0 C 
261 
262 C 
265 C 

264 r 

265 r 

266 c 

267 r 
266 C 

269 r 

270 r 

271 c 
2 7? r 
2 73 

274 

275 


X — abscissa array 
Y -- orhinaie array 
NLO — S1ARTII''IG INDEX OF PrGRrSSlOM 
^JHI -- EMDIMG index of REGRESSION 

slope: -- SLOPE OF y versus X 
YINT -- Y intercept 
ER -- ReiS error 

yOUT -- array OE SLOPE aND v INTERCEPT LINE 

( yout = slope»x + yint ) 

nIMENSION X(l)f Y(1),Y0UT(10 ?ii ) 

PAjA X0/-2./ 

XSDM=0. 

YS''N = P . 

XYSUf'lrO. 
y?.SUM=0 . 

calculate. SLOPE AND Y INTEKCEPf 
no 1 i=nlo,nui 

ySUM = XSUM + X( I )-Xl) 

YSLII5=YSUN + Y( I ) 

XYSUM=XYSUN+(X( I >-XO)*Y( X ) 
x2sUN=X?SUM-^ ( X n > -XO > **2 
FX = XSUM/(NhT-.NI.U+1 ) 
rY = YSUn/(Nf.I-NLO + i ) 

EXY = XYSUM/(Nhi-''*E041) 

EX2=X2SUfi/(NHI-NLO-*-l) 

SLOPE= (EXY-Ex*EY)/( Ex2-Ex*Fy > 

YIMT=EY-SL0PE*FX 

CALCULATE RMS ERROR 

FRRSUMsO, 

no 2 I=NLOfNHi 

youth )=SL0RE*(X< I )-X0)+YINT 
ERRSUrinERRSUM+I HOllT ( I ) - Y ( I ) ) *♦2 ) 
EK=ERRSUM/(Nhi-NLO+1) 

E R = SoRT(ER) 

RETURN 

ENO 


SUBROUTINE CLRHt YfNLO.NHlfSl.OPE,ER,YOuT ) 

CONSTRAINED LINEAR REGRESSION 

X -- ABSCISSA ARRAY 
Y — ORDINATE ARRAY 
MLO -- STARTING JNDEX 
NHI -- ENDING INDEX 
9L0PE — sSLOPE of Y VERSUS X 
ER RI^S fPROR 
YOUT -- array of SLOPE LINE 

PI MENS I ON X(1)«Y(1) fYOUT(in?4) 
rOMMON/SETl /NFRE.O»riHATE»TLF 
X(j=NLO»TLF 
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276 

277 
27« 

279 

2 R 0 

201 r 

202 
28.5 

204 

205 10 

206 

267 

268 
269 

290 r 

29J 

292 

295 

294 20 

295 

296 

297 

298 

299 r 

300 r 


X0=:ALOR10<Xtn 
YO = YOUT(NI.O) 

XSUI5=0. 

V2SUM=0, 
xY?un=o. 
calculate slope 
no in i=nld«nhi 

XSIJM = XSUM + X( I )-X0 
X2SUM=X2SUM+(X( I »"X0)**2 
XY9Un=XYSUM+ (X « I ) -XO ) ♦Y ( T ) 

pXsXSUM/tNHl-^JLO'fl ) 

E:X2 = X2SUfV ( NHl -Nl.0 + 1 ) 

FXY = XYSUP/(^HT-AJL0<-1 ) 
c;lOPE=:(EXY.£X*YO)/FX2 

calculate P^S FWPOR 
FHRSUM=n» 

no 20 I=(^LO,NhI 

YOUTI I )=SL0FE*( I ) -XO ) +yO 
EHRSUf1::EHPSUM+ ( TOUT ( 1 ) -Y (I ) ) **2 
fr=errsum/(Nhi-ailo+i) 

FR=SQRT (E R) 

return 

E I'jn 


72 



APPENDIX B 


REFRACTIVITY OF AIR IN THE PRESENCE OF 
WATER VAPOR AND OXYGEN 


Refractivity, N, can be expressed in general as 


N = 


S^Fl 


n;, 

w 


j(J 

1 


S-F‘i 

V 1 


N") 


ppm 


(25) 


where summations are over all spectral lines and 

N = frequency-independent refractivity (ppm) 

= line strength (KHz) , 

F = line shape factor (GHz“’^) 

N^ = non-resonant refractivity due to water vapor (ppm). 

1 q 

The frequency-independent refractivity is given by 


Ng = 77.67 + 71.7 =^5! 


p 

+ 3.744 X 10® (ppm) 


where P = dry air pressure (mbars) 

P = water vapor pressure (mbars) 
T^ = temperature {\) 


(26) 


Water Vapor . The resonance lines for water vapor are described 
hy 


p ' = i 
f 


fj -f 




1 




(GHz'^) 


(?7) 


pil _ f 
- j- 


Jl 


+ 




(GHz'^) 


(28) 


The only water vapor line relevant to this study is characterized 
by 


= 4.465 X 10^ (kHz) 
Y^ = 0.0861 (P + 4.8 P„)T‘°'® (GHz) 
fj = 22.235 (GHz) 


73 



The nonresonant contributions to refract ivity are 


w 


w 


= ?A32 P f 
w 

= l.T P P.,f 
w 



(ppm) 

m 

1. 3^-3. 5 

(ppm) 

f30) 


Oxyqen. Since there are numerous oxygen lines it is ^convenient 
to describe these lines by a single line characterized by'^ 


f -f 

r 

P V300 ^ 

2 

^2 ^ 

\ 1013. 25/V T } 


2 


/f9~f\2’ 


^2 


1 n 1 




. V Yg 7 . 


(. 

P 

^^oo 

1 V1013.257V T 7 

^2 

1 + 




W2 7 J 


(31) 


(32) 


along with the parameters 
= 5.22 kHz 
= 3.92 GHz 
= 60.0 GHz 

Derivatives of the Refractivity 

To compute the lower break frequency of the spectrum of am- 
plitude fluctuations it is desirable to obtain an expression for 
the derivatives of the real and imaginary parts of the refractivity 
with respect to temperature. 

a) The real part of the refractivity 
is found from Equation (25) as 


Re(N) = Nq + 


2 

I 

i=l 


S.Fl 
1 1 


+ N. 


w 


(33) 
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Each of the terms in the above equation is temperature dependent. 
The derivative of Re(N) with respect to temperature is found as 


3 Re(N) _ ^*^0 


3S 


+ F- 


1 


3F.' 3FA 3N' 

i . o C. , W 


3T ' '1 3T ^ 3T ^2 3T * 3T 

(34) 


where 


3N 

TT 


P P 

° = -77.67 ^ - 71.7 - 7.488 x 10^ ■— (ppm/k) 

T-%J 


3S 


3T 


- = -1.563 X 10^ P T"^*^ + ^^871 x 10^^ P 

w w 


(GHz/k) 




f 


3T 


f,+f 


^1-f 




j-y-^KGHz'^-k'’) 


9 Y 

gji = -0.05166 (P+4.8 PJT"^'® 


3F^ 

TT 


= -1.8x10 


5( 

2 

^2 


( P ) 

V1013.257 


1 + 


f -f 


(GHz'^-k‘^) 


3N;, 

W 

3T 


•2918 


(ppm/k) 


b) The imaginary part of the refract ivity is found from Equation 
(25) as 



The derivative of Im(N) with respect to temperature is found as 

.T i^\ as, 8FV 3 f;i 3N", 

^ =>^1 (ppm/k) (36) 


where 


fp" 

'^1 - ff 




(GHz‘’-k‘^) 


6 ( P ) 

-1.8x10^ V 101 3,26/ 


^ 2 "^ ^ "1 


^ = -3.85 PP (ppm/k) 
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APPENDIX C 


SPECTRAL SLOPE OF A RANDOM TELEGRAPH SIGNAL 


It is desirable to find the slope of the power density spec- 
trum of a random telegraph signal (which is similar to random tog- 
gling of the least significant bit of an A/D converter). 

From Papoulis , the autocorrelation function, R, of a random 
telegraph signal is exponentially damped as 

R{x) = (37) 


where t is the lag time and X is the occurrence rate. The power 
density spectrum, S(a)), is the Fourier Transform of R(t), 


S(o)l 


4X 

7 7 

4X‘^+(o- 


rs8) 


The spectral slope is 


6S{ixi) 

do3 


-8Xoj 

(4X^+0)^)^ 


f3«) 


An appropriate value of X can be obtained from Figure 28 of 
Chapter IV. There are approximately 110 bit changes in the 204.8 
second sample which gives a value of 110/204.8 = 0.54 for X, 


If the spectral slope is evaluated in the 0.25 Hz to 0.5 Hz 
region, X ^ w and the spectral shape is approximately 


dS(Ml (40) 

25u)^ 

Thus, the spectral slope of a waveform consisting of random 
bit toggling is -20 dB/decade, which is 6.7 dB/decade flatter than 
the spectral slope due to a Kolomogorov turbulence spectrum. 
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APPENDIX D 


SPECTRAL SLOPE OF STRONG FLUCTUATIONS 


Ishimaru^^ has considered the effect of strong fluctuations 
on the mutual coherence function, r. The mutual coherence func- 
tion is similar to the covariance function defined in Equation 
( 5) except that r is defined in terms of the specific intensity, 
I, instead of the log amplitude fluctuations, i.e.. 


r(x) = J ds (41) 

L 

where x is the spatial lag distance, s is the path vector, and 
the integration is along the path, L. 

If the turbulence is strong, x is large, and Taylor's hypoth- 
esis is valid then for a spherical wave 


rfx) 




(42) 


where 

F{x) = exp{-1.46k^(vT)®/^ / C^(x')dx'} , 

0 

2 

and V is the perpendicular wind speed, t is the lag time and C^ 
is the refractive index structure constant. 

Because Taylor's hypothesis is assumed, the product vt is 
slowly varying; furthermore, changes in vt do not have a signip'- 
cant effect on F(x). If it is assumed that vt < 10 then F(x) ^ 1, 
and the Fourier transform of r(x) is simply 

S((1)) = -TTO), 0)> 0 

which is the power density spectrum of the strong fluctuations. 

Thus, the slope of the strong fluctuations' spectrum is 0 dB/decade. 


78 



REFERENCES 


1. D. M. J. Davasirvatham and D. B. Hodge, "Amplitude Scintilla- 
tions on Earth-space Propagation Paths at 2 and 30 GHz," Re- 
port 4299-4, March 1977, The Ohio State University Electro- 
Science Laboratory, Department of Electrical Engineering; 
prepared under Contract NAS5-22575 for NASA-Goddard Space 
Flight Center. 

2. H. G. Booker and W. E. Gordon, "A Theory of Radio Scattering 
in the Troposphere," Proc. I.R.E., vol. 35, April, 1950, pp. 
401-412. 

3. F. VI liars and V. F. Weisskopf, "On the Scattering of Radio- 
waves by Turbulent Fluctuations of the Atmopshere," Proc. 
I.R.E.,'vol. ^3, October, 1955, pp. 123-1239. 

4. R. A. Silverman, "Fading of Radio Waves Scattered by Dielectric 
Turbulence," J, of Appl. Phys., vol. 28, April, 1957, pp. 
506-511. 

5. A. D. Wheel on, "Relation of Radio Measurements to the Spectrum 
of Tropospheric Dielectric Fluctuation," J. of Appl. Phys., 
vol. 28, June, 1957, pp. 684-693. 

6. Armel Picquenard, Radio Wave Propagation , John Wiley and Sons: 
New York and Toronto, 1974, Ch. 2. 

7. V. I. Tatarskii, Wave Propagation in a Turbulent Medium , R. 

A. Silverman, Transl., McGraw-Hill: New York, 1961. 

8. Akira Ishimaru, "Temporal Frequency Spectra of Multi -frequency 
Waves in Turbulent Atmosphere," I.E.E.E. AP-20, No. 1, January, 
1972, pp. 10-19. 

9. V. I. Tatarski, "The Effects of the Turbulent Atmosphere on 

Wave Propagation," transl. from Russian. Israel Program for 
Scientific Translations: Jerusalem, 1971. 

10. W. Lee Robert and Jeffrey C. Harp, "Weak Scattering in Random 
Media, with Applications to Remote Probing," Proc. I.E.E.E., 
vol. 57, No. 4, April, 1969, pp. 375-406. 

11. Akira Ishimaru, "Theory and Application of Wave Propagation 
and Scattering in Random Media," Proc. I.E.E.E., vol. 65, 

No. 7, July, 1977, pp. 1030-1061. 


79 



IP., 


P. A, Mandics, R. W. Lee, and A. T. Waterman, Jr., "Spectra 
of Short-term Fluctuations of Line-of-Sight Signals: Electro- 
magnetic and Acoustic," Radio Science, vol. 8, No. 3, March, 

1973, pp. 185-POl. 

13. Akira Ishimaru, "Fluctuations of a BeamWave Propagating through 
a Locally Homogeneous Medium," Radio Science, vol. 4, April, 

1969, pp. 255-305. 

14. R. H. Ott and M. C. Thompson, Jr., "Atmospheric Amplitude 
Spectra in an Absorption Region," I.E.E.E. Trans. Antennas 
Propagat. , vol. AP-26, No. 2, March, 1978, pp. 329-332. 

15. A. Papoulis, Probability, Random Variables, and Stochastic 
Processes , McGraw-Hill, New York, 1965, Ch. 9. 

16. R. B. Blackman and J. W. Tukey, The Measurement of Power Spectra , 
Dover Publications, Inc,, New York, 1958. 

17. Warren D. Brown, "A Model for the Refractive Index Structure 
Constant at Microwave Frequencies," Sandia Laboratories, 

New Mexico, February, 1977. 

18. W. I. Vogel, A. W. Straiton, and B. M. Fannin, "ATS-6 Ascending: 
Near Horizon Measurements over Water at 30 GHz," Radio Science, 
vol. 12, No. 5, September-October, 1977, pp. 757-765. 

19. H. J. Liebe and G. G. Gimmestad, "Calculation of Clear Air 
EHF Refractivitv," Radio Science, vol. 13, No. 2, March-April, 
1P78, pp. 245-251. 

PO. R. H. Ott and M. C. Thompson, Jr., "Characteristics of a Radio 
Link in the 55 to 65 GHz Range," I.E.E.E. Transactions on 
Ant. and Prop., AP-24, November, 1976, pp. 873-889. 

21. A. Papoulis, The Fourier Integral and its Applications , McGraw- 
Hill: New York, 1962, Ch. 12. 


80 



ACKNOWLEDGMENTS 


The authors wish to express special thanks to Dr. David M. 
Theobold for his numerous and fruitful ideas, especially those de- 
veloped in Appendices C and D, given to them throughout the course 
of this work. 

Thanks also go to Robert C. Taylor for his expert assistance 
during and after the experiment, to D. M. J. Devasirvatham for 
the preliminary data analysis which inspired this work, and to 
Dorothy M. McGinty for typing the manuscript. 


81 




